Electronic versus Phononic Friction of Xenon on Silver 
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Molecular dynamics simulations of a Xe monolayer sliding on Ag(OOl) and Ag(lll) are carried 
out in order to ascertain the microscopic origin of friction. For several values of the electronic contri- 
bution to the friction of individual Xe atoms, the intra-overlayer phonon dissipation is calculated as 
a function of the corrugation amplitude of the substrate potential, which is a pertinent parameter to 
consider. Within the accuracy of the numerical results and the uncertainty with which the values of 
the relevant parameters are known at present, we conclude that electronic and phononic dissipation 
channels are of similar importance. While phonon friction gives rise to the rapid variation with 
coverage, the electronic friction provides a roughly coverage-independent contribution to the overall 
sliding friction. 
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I. INTRODUCTION 

The microscopic origin of the sUding friction of thin 
films adsorbed on a metal surface has recently attracted 
considerable interest [MSI. The availability of refined ex- 
perimental methods makes it now feasible to in- 
vestigate the fundamental processes that contribute to 
the sliding friction on an atomic scale. On a metal sur- 
face, frictional energy is dissipated via two main channels: 
excitation of low-energy electron-hole pairs in the sub- 
strate |6|-|l2|] and excitation of phonons within the over- 
layer [n3-^|. Phonon excitation within the metal may 
also occur; however, this contribution usually has only a 
minor infiuence on the lateral sliding of adsorbed films. 

Detailed experimental information on the friction of Xe 
and Kr layers on noble metal surfaces was recently ob- 
tained by Krim et al. [^jH who performed quartz-crystal 
microbalance (QCM) measurements for a variety of cov- 
erages. In this technique, the frequency shift and broad- 
ening of the resonant line provide information on the 
number of adsorbed atoms and on their slippage dur- 
ing the lateral oscillation of the microbalance. The data 
were analyzed by various theoretical groups that, sur- 
prisingly, arrived at contradictory conclusions: For Kr 
on Au(lll), Cieplak et al. ||l|,|l4| argued that the ob- 
served friction is caused by lattice vibrations within the 
overlayer and that substrate-induced energy dissipation 
plays a negligible role. Dominant phonon damping was 
also inferred for Xe on Ag(lll) by Tomassone et al. [^ . 
On the other hand, Persson et al. argued that the 
phonon-associated friction for a full Xe monolayer on 
Ag(lOO) is small and that the measured friction is mainly 
of electronic origin. These opposite interpretations are 
puzzling since all three theoretical studies were based 
on molecular dynamics simulations for similar models of 
the rare gas/metal adsorption system. In addition, Lieb- 
sch [ [ll| , |r2| performed dynamical surface response calcu- 
lations within the time-dependent density functional ap- 
proach to determine the electronic friction of Xe atoms 



on Ag. The lateral component of the friction coefficient 
was shown to agree qualitatively with the experimental 
value for a full monolayer, suggesting that phonon pro- 
cesses play a minor role in the compressed phase. 

In view of these fundamental theoretical discrepancies, 
we decided to carry out molecular dynamics simulations 
for Xe on Ag, employing essentially the same overall 
model as previous authors |p^|-p^ . Key input parameters 
in these simulations are, in particular, the corrugation 
amplitude of the Xe/Ag interaction potential and the 
electronic friction of individual Xe atoms sliding parallel 
to the substrate surface. Since at present the corrugation 
amplitude is known from independent experimental data 
only with very appreciable uncertainty, a detailed under- 
standing of how its magnitude infiuences the final result 
is in our view crucial. In contrast to the earlier molecular 
dynamics studies, therefore, we do not choose a specific 
set of input model parameters but investigate the net 
sliding friction within a wide range of these parameters. 
In addition, we consider the dependence on the crystallo- 
graphic orientation by studying the friction of a Xe layer 
adsorbed on Ag(lll) and Ag(lOO), in order to check the 
importance of the lateral symmetry of the substrate po- 
tential. By carrying out such a systematic investigation 
we essentially cover the full range of overlayer/substrate 
models considered in the previous treatments. 

A striking feature of the Xe/Ag quartz-crystal mi- 
crobalance data by Krim et al. is the strong variation 
of the net sliding friction with Xe coverage. To facilitate 
the discussion, we show in Fig. 1 the measured slip time 
T as a function of Xe coverage Q on Ag(lll). The fric- 
tion coefficient fj is given by the inverse of the slip time. 
According to surface resistivity measurements for Xe on 
Ag the electronic friction coefficient depends weakly 
on coverage. It is evident, therefore, that the rapid cover- 
age dependence observed in the QCM data, in particular 
near completion of a monolayer, is caused primarily by 
phonons generated within the overlayer. This is not sur- 
prising since the probability of exciting phonons varies 
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greatly with inter-adatom spacing. With respect to this 
aspect of the data, there is httle dispute among the pre- 
vious simulations studies. 
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FIG. 1. Measured slip time r for Xe on Ag(lll), as a func- 
tion of coverage, after Ref. 4. In this paper we put forward 
arguments, based on extensive simulations, to show that the 
region roughly between the dotted lines can be understood 
as due to phonon originated friction, while above the upper 
dotted line a friction mechanism of electronic origin seems to 
prevail. 

The important question, rather, concerns the relative 
weights of electronic and phononic dissipation channels, 
especially at coverages near a full monolayer where the 
slip time reaches its maximum, i.e., where the sliding fric- 
tion of the incommensurate overlayer is smallest. Our 
results show that the effective friction parameter f] may 
be approximately written as 

fj = r]el+ Vph ■ (1-1) 

Here, ijei — V\\ is the parallel component of the 
single-atom electronic friction coefficient and the intra- 
adsorbate phonon contribution is roughly given by rjph = 
cuq, where uq is the corrugation amplitude of the Xe/Ag 
potential. A quadratic variation with uq agrees with the 
analytical and numerical results of Cieplak et al. . 
The coefficient c depends weakly on the crystal face and 
on the single-atom friction parameters r]\\ and r]±, but 
varies strongly with coverage. Evidently, according to 
the above expression, the measured sliding friction at 
any given coverage can be reproduced using several com- 
binations of the parameters uq and r]ei. For instance, by 
reducing the single-site electronic friction ?7|| , one is able 
to match the data if the substrate potential is assumed 
to be sufficiently strongly corrugated. Conversely, if the 
electronic friction is chosen large enough, little additional 
phonon damping is required to reproduce the same net 
fj. Constraints on the range of acceptable parameter val- 
ues are implied by the observed variation of the slip time 
with coverage, in particular, by the steep slope near com- 
pletion of the first monolayer. 



The molecular dynamics simulations for Xe on Ag, 
which we describe below, suggest that the characteris- 
tic coverage dependence of the observed slip time and its 
overall magnitude cannot be understood only in terms of 
phonon dissipation. Although there is considerable un- 
certainty both in the calculated and measured friction 
coefficients, the most likely scenario is the presence of 
a roughly coverage-independent electronic contribution 
of about T « 2 ns (ry|| sa 0.5 ns~^) and an additional, 
strongly coverage-dependent phonon friction, which can 
reduce the slip time to about r « 1 ns (the net slid- 
ing friction then increases approximately to 77 « 1 ns""'^). 
Thus, roughly speaking, we associate the slip time in the 
region between the horizontal dotted lines in Fig. 1 with 
phonon processes, while the portion above the upper line 
is related to electronic energy dissipation. Remarkably, 
the electronic friction deduced from this analysis is in 
qualitative agreement with independent surface resistiv- 
ity data for Xe on Ag ||] and with theoretical estimates 
of this dissipation channel pd| , p^ . Thus, the friction of a 
full monolayer is essentially of electronic origin, while the 
rapid increase of the slip time before reaching the com- 
plete monolayer is due to excitation of phonons within 
the Xe overlayer. 

This paper is organized as follows: Section II summa- 
rizes the theoretical results available for the electronic 
friction of Xe atoms on a metal surface. Section III 
presents the model employed in the molecular dynam- 
ics simulations. In particular, we specify the interaction 
potentials and the friction felt by a single Xe atom. Var- 
ious technical aspects of the calculation of the phonon 
friction are also discussed. Section IV contains the re- 
sults, with special emphasis on the variation of the net 
sliding friction with the corrugation of the substrate po- 
tential and its dependence on the electronic friction. Our 
conclusions are given in Section V. 



II. ELECTRONIC FRICTION 

We consider first the electronic contributions to the 
sliding friction of individual Xe atoms above a metal sur- 
face. In a physisorption system, these may arise due to 
the van der Waals attraction and the Pauli repulsion. In 
the case of Xe, some dissipation might also be associated 
with the formation of a weak covalent bond. 

It is well known that the instantaneous mutual po- 
larization between neutral species gives rise to the long- 
range van der Waals attraction. If an atom moves relative 
to a metal, the induced surface charge lags behind, caus- 
ing kinetic energy to be dissipated. Thus, as pointed out 
long ago by various authors , the excitation of low- 
frequency electron-hole pairs in the surface region con- 
tributes to the friction of neutral atoms. If the adatom 
velocity is small compared to the Fermi velocity of the 
metal electrons, the friction coefficient may be expressed 
as@,| 
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V^ = -— lim - [d^r [(fr'VJr)lmx{r,r',uj)Vi(r'). 
m i^-^o uj J J 

(2.1) 

The index i = x, z refers to parallel or perpendicular mo- 
tion of the atom, x is the many-body response function 
of the metal, and m is the adatom mass. The effective 
interaction Vi (r) is approximately given by 



r-d| 



(2.2) 



Here, a(0) is the static polarizability of the adatom and 
d = (0, 0, d) is the location relative to the jellium edge. 
Since Vi (r) does not satisfy the Laplace equation, it can 
only be approximately represented in terms of a super- 
position of evanescent plane waves of the form p(] 



■rii+gz 



(2.3) 



With this expansion we can write r/i as 

V^ = -l^ l™n 1 1^5(9,'^) , (2.4) 

m Ztt ij^o UJ ^ — ' " 

qii 

where the surface response function g(q,Lu) is defined as 

(2.5) 



gil,^) = / dz e*^ ni{z,q,Lj) 



and ni{z,q,uj) is the electronic surface density induced 
by an external potential given by — — e"?i ''11+'?^. 

In the case of Xe atoms physisorbed on Ag(lll) at a 
temperature T = 77.4 K, the distance d = 2.4 A. The 
static polarizability of Xe is a(0) = 4.0 A'^. Using the 
density functional results for ni{z,q,uj) and g(q, uj), we 
find - 0.34 ns-i and rij^ - 1.74 ns'^ Q^. 

We note that the theoretical value of ?7|| given above 
presumably represents a slight overestimate since the Xe 
atom is not completely outside the range of th e electronic 
density profile as implicitely assumed in (2.4). Moreover, 
as a result of the hybridization between the Ag s and 
d electrons, the surface polarizability of real Ag may be 
slightly smaller than that of the corresponding jellium 
model [|l^. We also point out that, according to (2.2), 
the friction coefRcient varies asymptotically like iJdP, 
i.e., it is rather sensitive to the precise location of the 
adatom above the metal surface. 

The friction coefficient arising from the van der Waals 
attraction derived above is significantly larger than the 
contribution due to the Pauli repulsion which was esti- 
mated by Persson [|j to be about ri\\ ~ 0.06 ns~^. In 
addition, for Xe there may exist some friction due to 
chemical effects resulting from the broadening of the Xe 
6s level. Estimates of this mechanism yield ^ ?7|| ~ 
0.15 ns"^ 



On the basis of these theoretical calculations and esti- 
mates, we conclude that the total electronic contribution 
to the lateral friction of single Xe atoms on Ag is ap- 
proximately 0.5 . . . 0.6 ns""'^. At full monolayer coverage, 
this friction is presumably slightly smaller because of the 
weakening of the Xe/Ag bonds. An electronic contribu- 
tion of this magnitude is in excellent agreement with the 
surface resistivity measurements for Xe on Ag 

According to our molecular dynamics simulations at 
T = 77.4 K, the distance of Xe atoms normal to the sur- 
face varies about 0.1 A around the equilibrium distance 
d — 2.4 A. For the asymptotic van der Waals attrac- 
tion, such a variation implies a variation of 77|| roughly 
between 0.2 and 0.5 ns^^. For the Pauli and covalent- 
bond contributions the variation is exponential due to 
wave function overlap, with an exponent determined by 
the work function. This leads to a variation of r;|| of 
about 20 %. Thus, for some Xe atoms, the electronic 
friction is larger, for others smaller than the equilibrium 
value 0.5 . . . 0.6 ns~^ quoted above. Because of this un- 
certainty the simulations are carried out for a sufficiently 
wide range of electronic friction coefficients. 



III. MODEL FOR PHONON FRICTION 

In this section we focus on the evaluation of the friction 
induced by the excitation of phonons within the adsorbed 
layer of Xe atoms. The single-atom frictional proper- 
ties are assumed to be known and are used as inputs in 
the molecular dynamics simulation. We specify first the 
intra-adsorbate and adsorbate-substrate potentials and 
then discuss details concerning the molecular dynamics 
simulations. 



A. Potentials 



As in previous work |13|-|16[|, the interaction between 
Xe atoms is expressed as a sum of Lennard- Jones pair 
potentials 



v{r) 



(3.1) 



where e = 19 meV is the well depth and rp = 4.54 A is 
the inter-particle spacing at the potential minimum. The 
total Xe-Xe interaction potential then takes the form 



(3.2) 



The adsorbate-substrate interaction potential is assumed 
to be given by 



U{r) = Eo 



uoe 



-a(z-zo) 



'u{x,y) 



(3.3) 
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(3.4) 



where for Ag(lOO) 

u{x, y) = 2 + cos{kx) + cos(fcy) 
whereas, in the case of Ag(lll), 



u{x, y) = 1.5 + cos ^2fcy/V3^ + cos {k{x + y/ %/3) 

+ cos(fc(x-y/V3)) . (3.5) 



Here, k = 27r/a and a = 2.89 A is the spacing be- 
tween neighboring Ag atoms. The unit vectors of 
the substrate lattice are (a,0), (0, a) for Ag(lOO) and 
(a, 0), (a/2, ay/2,/2) for Ag(lll). The binding energy for 
Xe on Ag is Eo = 0.23 eV. The decay constant a may be 
deduced from the frequency of the perpendicular vibra- 
tion ||l5|, i.e., a = 0.72 A~^. The decay constant a' of 
the Fourier component should in principle be larger than 
a. However, this difference is not very important and we 
choose a' = a. The corrugation factor u{x, y) is assumed 
to vanish at hollow adsorption sites. The potential bar- 
rier between neighboring hollow sites is 2uq for the (100) 
face and 0.5uo for the (111) face. The total barrier be- 
tween hollow sites across the substrate lattice sites is 4wo 
for the (100) face and 4.5wo for (111). 

The adsorbate-substrate interaction specified in 



Eq. (3^) has the same form as the one used by Persson 
and Nitzan |l^. Cieplak et al. |13 1^ and Tomassone et 
al. | [T^ , on the other hand, used the expression derived 
by Steele ||l8| for a sum of pair potentials. Since this 
potential is much more corrugated than the true inter- 
action between rare gas atoms and metal surfaces, the 
coefficient of the non-vanishing Fourier component was 
drastically scaled down. The functional form of this po- 
tential also differs from expression ( |3.3D . However, the 
precise shape near the potential minimum should not be 
crucial, as long as the overall corrugation is the same. 

Since the parameter uq is decisive for the excitation of 
phonons within the overlayer, we do not choose a particu- 
lar value. Instead, we calculate the net sliding friction as 
a function of uq in order to illustrate its sensitivity to the 
adsorbate-substrate potential. In this manner, we cover 
the various substrate potentials considered previously. 



B. Simulations 

There are several methods of extracting the friction 
parameter from the molecular dynamics simulations, 
(i) One can apply a constant lateral force F to all the 
adatoms, while keeping the substrate atoms at rest. Af- 
ter an appropriate number of time steps, the steady-state 
velocity v of the overlayer is determined. If the friction 
is viscous, V satisfies the linear relation 



mrjv 



(3.6) 



(ii) To simulate the quartz-crystal microbalance measure- 
ment, one may let the substrate atoms oscillate later- 
ally and determine the induced lateral vibration of the 
adatoms. This method was also used by Cieplak et 
al. ]l^ , |l^ ; (iii) One can apply a constant force to the 
adatoms up to a specific point in time and then switch it 
off. If the friction is viscous, the decay of the steady-state 
velocity of the adsorbate is proportional to e"***; (iv) Fi- 
nally, one may determine fj from the thermal equilibrium 
autocorrelation function of the center of mass velocity of 
the overlayer. The latter two methods were employed by 
Tomassone et al. [ p^ . 

In the present work, we also apply a constant lateral 
force F to the adatoms and determine the steady-state 
velocity that is reached after a sufficiently long time. The 
simulations are based on the Langevin equation 



jnrjri 



dU 
dr. 



dV 
dr,. 



(3.7) 



where fi is a stochastically fluctuating force describing 
the effect of the irregular thermal motion of the substrate 
on the i*'' adatom. The components of are related to 
the microscopic friction r/ via the fluctuation-dissipation 
theorem 



(3.8) 



where m is the adatom mass. This procedure was used 
by Cieplak et al. ^3 14 1 and by Persson and Nitzan ; 



The microscopic friction tensor r]ap is assumed to be di- 
agonal, with independent elements r/|| and rj±. 

To solve Eq. (3.7) the time variable is discretized with 
a step duration 0.01 ^Qj where the natural time unit is 
^0 — ^0 (m/e)^/^ = 3.84 ps. The simulations are carried 
out using the integration procedure suggested by TuUy 
et al. p9| |. Typically, thermalization was reached after 
10^ — l(rto. The drift friction fj was derived using the 
definition fj — F/{m{v)), where the external force F was 
applied in the z-direction on each adsorbate, resulting in 
the (time averaged) drift velocity (u). On Ag(lOO), the 
basic unit was taken as a square containing 12 x 12 or 
24 X 24 substrate atoms. Similar unit cells were employed 
for Ag(lll). These cells were then repeated assuming pe- 
riodic boundary conditions. 

Except for the form of the substrate potential, the 
model outlined above is very similar to the ones used in 
earlier simulations 13-1^1. The main difference between 
the treatments Of Refs. |l^-jlq| resides in the interpre- 
tation of the microscopic friction of a single Xe adatom. 
In the work by Cieplak et al. this friction plays 

the role of a thermostat that allows to establish a con- 
stant temperature T within the overlayer. Since the di- 
rection of this single-site friction is chosen orthogonal to 
the direction of the external force, it was claimed not to 
affect the net sliding friction of the overlayer. Instead, 
in the work by Persson and Nitzan |jl^ this thermostat 
was given a specific microscopic origin. In particular, t^h 
coincides with the lateral friction of individual Xe atoms 
due to excitation of electron-hole pairs in the substrate, 
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i.e., r^ii = rjei- According to the theoretical estimates and 
the resistivity measurements discussed in Section 11, this 
coefficient should be about r/|| w 0.5 ns^^. Finally, in 
the work by Tomassone et al. ||l^ , the simulations were 
carried out in the absence of a thermostat or any other 
single-site friction forces. 

In principle, electronic processes also contribute to the 
friction of the perpendicular motion of single adatoms. 
As pointed out in Section 11, density functional calcu- 
lations yield 77^ ~ 1.74 ns^^. However, this contribu- 
tion is much smaller than the normal friction induced by 
phonon excitation, which was estimated by Persson and 
Nitzan to be roughly ryj^ = 260 ns~^. Because of the 
present uncertainties in the evaluation of both and 
77j_, we have carried out simulations for a whole range 
of these parameters, in order to illustrate the sensitivity 
of the net sliding friction to the microscopic processes 
governing the behavior of individual Xe atoms. 



IV. RESULTS AND DISCUSSION 



In this section we first discuss the sliding friction of 
a monolayer of Xe atoms on Ag and then address the 
coverage dependence. In the final part, we compare 
these results with the QCM data. The temperature of 
T = 77.4 K, which we adopt, corresponds to the one used 
in the measurements by Krim et al. We define the 
coverage = 1 as the fully compressed monolayer pc| ] 
with density Ua = 0.0597 A'^. On Ag(lOO) this amounts 
to TV = 72 Xe atoms within the basic 12a x 12a unit cell 
(N = 288 within the 24a x 24a cell). The uncompressed 
monolayer has a slightly lower density, Ua = 0.0565 A~^, 
i.e., e = 0.94 {N = 68 atoms within the 12a x 12a unit 
cell). The remaining parameters required in the simula- 
tions, which we report below, have been specified in the 
preceding sections. 



A. Corrugation dependence of the sliding friction for 
an uncompressed Xe layer 



In Fig. 2 we display typical results of the center of mass 
velocity, in the direction of the applied force, as a func- 
tion of time. The examples are for = 68 on the (100) 
surface and the constant external force on each of the ad- 
sorbates F — 0.001 e/rp is apphed after 100 <o as indi- 
cated by the arrow. The single-atom friction coefficients 
are t^h — and 1.3 ns"'^, with ri± — 260 ns~^. The 
corrugation amphtude is uq — 0.95 meV and 1.9 meV. 
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FIG. 2. Center of mass velocity (in units A/ps) of a 
Xe layer on Ag(lOO) as a function of time. The coverage 
corresponds to an uncompressed monolayer at Q = 0.94 
or Ua = 0.0565 atoms/A'^. The arrow indicates the 
time {t — 100 to — 384 ps) at which the external force 
F = 0.001 e/ro = 4.2 x 10"^ meV/A is applied. Upper curves 



are for T^y — 0, lower curves correspond to 
(a) uo = 0.95 meV; (b) uq = 1.9 meV. 



^11 



1.3 ns 



One of the conclusions of the present work can be de- 
rived directly from this figure. In Fig. 2a, for a relatively 
low corrugation amplitude uq = 0.95 meV, the effect of 
the single-site friction coefficient ryu is evident: suppress- 
ing it, the center of mass velocity increases and so do the 
fluctuations. But the effect of the suppression is attenu- 
ated when the corrugation amplitude is bigger, as is the 
case in Fig. 2b for uq = 1.9 meV. Thus, the relative im- 
portance of phononic or electronic dissipation channels 
depends crucially on the magnitudes of uq and rj^^ . 

The validity of the linear response regime can be 
checked in Fig. 3, where we show the center of mass 
velocity {v) , time averaged over 4000 to after the tran- 
sient period following the application of the external force 
F, as a function of F. The approximate linear relation 
F — mf]{v) is apparent in all cases, indicating that the 
system indeed obeys a viscous force law. The full lines 
represent weighted linear regressions of the (v) ^ F re- 
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lation, which we use to derive the net sliding friction rj. 
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FIG. 3. Average center of mass velocity (v) (in units of 
ro/to) as a function of external force F (in units of e/ro) 
showing the linear viscous regime for different values of the 
input parameters t^h and uq. Each point was obtained by time 
averaging over a long steady state run and over two indepen- 
dent runs. Error bars were estimated from the fluctuations of 
the center of mass velocity as can be seen in Fig. 2. Straight 
lines are linear fits of points accounting for error bars. The 
values of -q are as follows (from top to bottom): 0.0023/to, 
0.0045/to, and 0.010/to. 

Fig. 4 summarizes our results for as a function of 
the corrugation amplitude uq and for different values of 
the lateral electronic single-atom friction r]\\ . The overall 
sliding friction is approximately of the form 



r] ^ 1]el + CUq , 



(4.1) 



with riei 
tion r]ph 



T^ll . The intra-adsorbate phonon contribu- 
cuq varies quadratically with the corruga- 
tion amplitude uq and the coefficient c depends weakly 
on the coefficient ?7|| and on the crystal structure of the 
substrate. This can be seen by comparing the results 
shown in Fig. 4a and 4b for the Ag(lOO) and Ag(lll) 
surfaces, respectively. Although the overall variation of 
fj is similar for both faces, the coefficient c is slightly 
larger for Ag(lll) than for Ag(lOO), indicating that fj is 
governed by the total barrier height (4.5 mq for (111) vs. 
4 Mo for (100)), rather than by the corrugation barrier be- 
tween neighboring hollow sites (0.5 uq for (111) vs. 2uq 
for (100); see previous section). This trend is plausible 
since, at monolayer coverage, the Xe atoms are not free 
to move along paths between nearest potential minima 
but are forced across potential maxima. We note here 
that, in the ?7|| — quadratic variation of rj with 

Mo was also found by Cieplak et al. [^,0 . 
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FIG. 4. Net sliding friction of Xe layer on Ag as a function 
of corrugation amplitude uq. The coverage is O = 0.94. The 
results are shown for three values of the lateral single-atom 
electronic friction, indicated by the symbols: r^n = rjsi = 0, 
0.0025 /to, and 0.0050 /to. tjx = 1.0 /to in all cases, (a) (100) 
surface; (b) (111) surface. 



In order to demonstrate the weak dependence of the co- 
efficient c on r/|| , we collapse in Fig. 5 the data obtained 
for different ry|| into one curve by plotting 



Vph = V~V\\ 



(4.2) 
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The full lines correspond to a quadratic fit of the data 
points. According to Figs. 5a and 5b, we deduce for the 



(100) and (111) faces c« O.Gip 
and c « 0.8i^^e-2 ~ 0.56 ns-imcV" 



0.42 ns^imeV 
, respectively. 
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FIG. 5. Phonon contribution to the sliding friction of Xe 
layer on Ag obtained as — r\\^ , as a function of corrugation 
amplitude uq. Symbols as in Fig. 4. Full lines are quadratic 
fits accounting for errors of points. Error bars are the same as 
Fig. 4 but are suppressed here for clarity, (a) (100) surface; 
(b) (111) surface. 

The results shown in Figs. 4 and 5 correspond to 
= l/<o = 260 ns^^. Simulations for different val- 



ues of Tyj^ in the range 26 ns""'^ < ?7± < 260 ns""'^ yield 
nearly the same overall sliding friction, confirming the 
results of previous studies which suggested that the fric- 
tion in the direction of the external force is essentially 
independent of the thermostat or the single-site friction 
in the orthogonal direction |p^|-|l5| . 

The simulations discussed so far were carried out for 
68 Xe atoms within substrate unit cells corresponding to 
12a X 12a for Ag(lOO) and 12a x 14a^/3/2 for Ag(lll). 
Remarkably, we observed pronounced finite-size effects 
associated with the small dimension of these unit cells. 
This is illustrated for Ag(lOO) in Fig. 6, where r]^h is plot- 
ted as a function of uq. In general the Xe atoms form ap- 
proximate hexagonal layers. Their orientation, however, 
tends to prefer two angles: with rows of Xe atoms nearly 
parallel to the x-axis, or nearly parallel to the y-axis. We 
denote these phases by a and /3, respectively. Since these 
orientations persist even in the absence of the substrate 
corrugation (mq = 0) and for vanishing external force, it 
is clear that they are a consequence of the limited size 
of the unit cell. The reason is that it is not possible to 
accomodate truly incommensurate hexagonal overlayers, 
at arbitrary angular orientations, in a rectangular unit 
cell that is periodically repeated. 
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Mq (meV) 

FIG. 6. Phonon contribution to the sliding friction of Xe 
layer on Ag(lOO) obtained as 77 — t^h, as a function of corru- 
gation amplitude uq. Symbols and quadratic fits as in Fig. 5. 
Full symbols and full line: a-phase, open symbols and dashed 
line: /3-phase, corresponding to different angular orientation 
of Xe monolayer. 

In order to reduce these finite-size effects, we have per- 
formed simulations for much larger substrate cells (e.g., 
288 Xe atoms in a 24a x 24a cell on Ag(lOO)). Quali- 
tatively, the net sliding friction for these systems agrees 
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rather well with the earlier results, but the tendency for 
orientational alignment of the Xe layer is greatly reduced. 



B. Coverage dependence of sliding friction 



As mentio ned above, the coefBcient c of the quadratic 
term in Eq. (4.1) depends strongly on overlayer coverage. 
We illustrate this point in Fig. 7, where the net slid- 
ing friction of the uncompressed monolayer (0 = 0.94) 
is compared with that for a fully compressed Xe layer 
(9 = 1.0) and for a less dense layer (6 = 0.85). 
The compression is seen to lead to a greatly reduced 
phonon-induced contribution to fj. On Ag(lOO) as well 
as Ag(lll), the magnitude of c is reduced by more than 
a factor of two. A strong reduction of dissipation via 
phonons upon compression is to be expected since the Xe 
atoms are much less free to vibrate about the equilibrium 
positions. The opposite trend is found if the coverage is 
reduced below that of the uncompressed monolayer. For 
Xe on Ag(lOO) at coverage 9 = 0.85, the coefficient c is 
about 1.25 ns~^meV^^, i.e., three times larger than for 
9 = 0.94. An even more dramatic increase of c is found 
on Ag(lll). The strong enhancement of fj is caused by 
the higher probability of exciting atomic vibrations in the 
less dense overlayer phase. 
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FIG. 7. Net sliding friction r\ of Xe layers on Ag 
as a function of corrugation amplitude mq for coverages 
e = 0.85, 0.94, 1.0. (a) (100) surface; (b) (111) surface. 
r\\\ — 0.65 ns~^. 

Concerning the overall variation of r\ph hi the coverage 
range 9 = 0.85 . . . 1.0, on Ag(lOO) for = 0.95 meV 
the variation is about 1 ns~^, while for uq — 0.5 meV 
this variation diminishes to less than 0.3 ns~^, and for 
1*0 — 1.4 meV it is larger than 2 ns^^. On Ag(lll), in 
the same coverage range, uo = 0.3 meV yields a vari- 
ation of 77p^ of about 1 ns^^, whereas uq — 0.2 meV 
(uo — 0.4 meV) lead to much smaller (larger) variations. 
Therefore ryp^ is sensitive to the substrate topology. 

In view of the great sensitivity of ?7ph(9) with respect 
to the corrugation amplitude, it would be of interest 
to have more accurate independent experimental and/or 
theoretical determinations of the substrate/adsorbate po- 
tential. In principle, the corrugation amplitude may 
be estimated from the lateral vibrational frequency t^n 
of adsorbed Xe, which can be detected using inelastic 
He scattering. For Xe on Cu(lll), one finds ~ 
3 ± 1 cm~^ Using the relation uq = muy^jk'^^ 

one obtains values of uq in the ap prox imate range be- 
tween 0.5 and 2 meV. From Eq. (4.1) it follows that 
r\ph varies with the fourth power of This is born 
out by our simulations which show that such barriers 
yield a phonon-induced friction in a very wide range, 
•r]ph ~ 0.1 .. . 5.0 ns^^, if the coverage is varied between 
1.0 and 0.85. 



C. Comparison with Experiment 

As pointed out in Section I, the quartz-crystal mi- 
crobalance measurements by Krim tt al. M for thin Xe 
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films on Ag(lll) show that the sHp time r = l/fj ex- 
hibits a characteristic dependence on Xe coverage (see 
Fig. 1). Close to one monolayer, the slip time reaches 
a minimum of about t « 0.8 ns and then increases 
to about 2.0 ns within a coverage range of less than 
0.3 monolayers. In this section we focus on this cov- 
erage range (roughly the "compression region") since it 
provides the most challenging aspect of the data. The 
interpretation of the slip time at submonolayer coverages 
in terms of simulations based on small repeated cells is 
less reliable because they do not account for the forma- 
tion of islands. Also, the presence of defects and steps 
which are not included in the theory then plays a larger 
role. On the other hand, beyond full monolayer coverage, 
the observed slip time shows a much weaker coverage de- 
pendence related to additional phonon dissipation in the 
partly filled second monolayer. 

In Ref. it is stated that the minimum of t 

corresponds to the uncompressed monolayer at tIq « 
0.0563 and that the subsequent steep rise reflects 
the suppression of intra-adsorbate phonon excitations as 
the Xe layer becomes more compressed (as noted above, 
the compressed phase is only 6 % denser than the un- 



compressed layer |20 .) This coverage assignment implies 
that the slip time increases (i) when the coverage is re- 
duced below that of the uncompressed monolayer; and 
(ii) when the coverage is increased beyond that of the 
compressed layer. However, both of these implications 
are implausible since in general incommensurate solid 
layers slide more easily than fluids ||l^-[l5[|. Thus, the 
slip time should decrease rather than increase (i) when 
the coverage is lowered below that of the uncompressed 
layer and (ii) when it is increased beyond that of the com- 
pressed layer, i.e., when the second layer is beginning to 
adsorb. Since absolute coverage calibration in the ex- 
periment is difficult, there exists a clearly appreciable 
uncertainty concerning the measured sliding friction of a 
monolayer. In fact, it would appear more plausible to 
associate the minimum of the slip time not with the un- 
compressed layer but with a slightly lower coverage. Part 
of the observed steep rise of t would then correspond to 
the reduction of intra-overlayer phonon processes upon 
reaching the uncompressed layer coverage, and the re- 
maining rise should reflect the additional phonon reduc- 
tion due to compression. The maximum of the observed 
slip time should approximately correspond to the fully 
compressed phase. According to these arguments, the 
uncompressed monolayer might have r « 1.3 ... 1.6 ns, 
giving ^ « 0.6 . . . 0.8 ns"""^. 

In Fig. 8 the phonon friction ripfi{Q) is plotted as a 
function of coverage, in the range O = 0.85 ... 1.0, for 
several values of the corrugation amplitude uq. For ease 
of comparison with the experimental results, we include 
in Fig. 8(b) the measured 77 = 1/r as a function of cov- 
erage in the region near one monolayer where r exhibits 
the characteristic positive slope. To account for the un- 
certainty of the measured coverage, and for the expected 
behavior of the sliding friction as discussed above, a shift 



of the data to lower coverages, by about 0.1 . . 
layer, should be allowed for. 
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FIG. 8. Calculated phonon friction rjphiQ) of Xe on Ag as 
a function of coverage Q for several corrugation amplitudes 
Uq. (a) (100) surface; (b) (111) surface. The dashed lines 
are guides to the eye. The error bars indicate the uncertainty 
due to fluctuations. Panel (b) also shows the measured slid- 
ing friction 77 = 1/r in the compression region where the slip 
time T exhibits the characteristic positive slope (see Fig. 1). 

The coverage dependence of fj in the range shown can 
be characterized by two important features: the overall 
magnitude and the slope. The comparison of fj with the 
calculated results shows that none of the functions riph{Q) 
obtained for different corrugation amplitudes uq pro- 
vides an adequate description of the data. However, for 
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uq = 0.19 meV, the calculated r]ph has roughly the right 
slope even though its magnitude at full monolayer cover- 
age is far too small. If uq is increased to 0.3 . . . 0.4 meV 
the slope of Tjph{&) becomes rapidly larger than that of 
the measured f] while the magnitude near monolayer cov- 
erage is still too small. Only when uq is increased to much 
larger values does rjph, in the compressed phase, match 
the observed fy, but the slope of r]ph{&) then is much 
larger than observed. 

As shown in Fig. 8(a), a similar picture is found on 
Ag(lOO). Although the detailed dependence of riph{Q) 
on the corrugation differs appreciably from the one ob- 
tained for the (111) face, near monolayer coverage, rjph 
remains less than 0.1 ns~^ as long as mq < 1 meV, with 
a moderate slope for 8 < 1. Only if uo is increased to 
2 meV we find rjpfi « 0.5 ns^^ in the compressed phase. 
The slope of r]ph{Q), however, then quickly becomes ex- 
tremely large. 

Thus, our results imply that it is not possible to simul- 
taneously describe the magnitude and slope of the ob- 
served variation of fj^Q) exclusively in terms of phonon 
dissipation. Only by adding an electronic friction con- 
tribution (which in the narrow range shown in Fig. 8 
can be assumed to be independent of coverage) does 
the sum rjei + rjphiQ) reproduce the measured variation 
of ?7(0). The comparison with the data suggests that 
T]ei « 0.5 ns^^ together with uq « 0.2 meV (total barrier 
height « 0.9 meV) provides a qualitative description of 
the sliding friction measured on Ag(lll) — if one allows 
a shift of the data of 0.14 monolayer to lower coverages, 
of course. As discussed in Section II, an overall single- 
atom electronic friction coefficient of this magnitude is 
in agreement with surface resistivity measurements and 
independent theoretical estimates. 

Interestingly, our findings for Xe/Ag differ from those 
on Kr/Au. As recently shown by Robbins and Krim 
p2| , the Kr/Au QCM data up to monolayer coverage 
can be fitted rather well assuming negligible electronic 
friction. A finite electronic friction would require a cor- 
respondingly smaller substrate corrugation uq which, in 
turn, would yield too low values of the phonon friction 
at submonolayer coverages. We note, however, that the 
analysis of the low-coverage data is uncertain since the 
simulations are limited to repeated finite cells and do not 
include the effect of defects or steps, nor do they properly 
describe the formation of islands which most likely occur 
under experimental conditions. Thus, if one focusses on 
the slip time in the less problematic region near mono- 
layer coverage, i.e., roughly 0.85 < 6 < 1, the analysis 
in Ref. indicates that the inclusion of a finite elec- 
tronic friction and a slightly reduced corrugation uq is 
not incompatible with the data. 

In addition, we point out that electronic friction of Kr 
on Au indeed ought to be weaker than for Xe on Ag for 
the following reasons. Since the Au d bands lie only 2 eV 
below the Fermi level, in contrast to 4 eV for Ag, the 
larger s — d hybridization diminishes the polarizability 
of the s electron tails decaying into the vacuum region 



and thereby reduces the probability of exciting substrate 
electron-hole pairs via adsorbed particles. The larger 
work function of Au (5.3 eV compared to 4.8 eV) also 
contributes to a stiffer, less polarizable electronic density 
tail JlT| , p^ . Moreover, since the atomic polarizability of 
Kr is about 40 % smaller than that of Xe ||2^] , the leading 
van der Waals part of the friction coefficient should be 
much smaller than for Xe on Ag. For the same physi- 
cal reasons, the contributions associated with the Pauli 
repulsion and the covalency of the adatom bond should 
also be smaller. The overall electronic coupling of Kr 
on Au should therefore be significantly weaker than for 
Xe on Ag. Thus, great caution is required when compar- 
ing electronic excitation processes on different adsorption 
systems — even in the case of rare gas atoms on noble 
metal surfaces. 

The comparison of Figs. 8(a) and (b) suggests that, 
for the same corrugation amplitude uq, Ag(lll) yields a 
larger phonon friction than Ag(lOO). This trend seems at 
first surprising since the close-packed (111) face is gen- 
erally regarded as the smoothest of the low-index faces 
of an fee crystal. On the other hand, as pointed out 
in Section III, the potential barrier between neighboring 
hollow sites is 2uo for the (100) face and 0.5uo for the 
(111) face, while the total barrier between hollow sites 
across the substrate lattice sites is 4uo for the (100) face 
and 4.5mo for (111). Thus, the (111) surface has sharper 
potential maxima than the (100) face. Since atoms of 
an incommensurate overlayer sample not only the poten- 
tial minima but the entire substrate surface, it is indeed 
plausible that Ag(lll) causes more Xe intra layer phonon 
friction than Ag(lOO). According to this, the topology of 
the substrate is a relevant feature if one want to make a 
comparison with the experiments. 

In their analysis of an uncompressed Xe layer on 
Ag(lOO), Persson and Nitzan ||l^ used uq = 0.95 mcV 
and obtained ?/ph ~ 0.01 ns^^ (for 7y|| = 0.62 ns~^, they 
find ?/||/?7 = 0.98±0.04). From this result, they concluded 
that the net friction at monolayer coverage is mainly of 
electronic origin. Instead, for the same parameters, we 
find riph ~ 0.3 ns^^, suggesting a phonon friction roughly 
half as large as the electronic contribution. 



On the other hand, Tomassone et al. |16| assumed a 
total barrier height of 2 meV on Ag(lll), implying in 
our notation uq = 0.45 meV. On the basis of the simu- 
lations it was concluded that phonon friction dominates 
over possible electronic contributions. For an uncom- 
pressed Xe layer, these authors obtained a phonon fric- 
tion rjph ~ 1.8 ns""'^, which is considerably larger than 
what we find for the same parameters. According to 
Fig. 2 of Ref. , in the compression region, the calcu- 
lated slip time for uq = 0.45 meV increases extremely 
rapidly in agreement with our results. This increase, 
however, is much steeper than what is seen in the ex- 
periment. Thus, a weaker corrugation, combined with a 
finite electronic friction might, in fact, provide a better 
fit of the data in this coverage range. The analysis of 
the submonolayer data is again more complicated for the 
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reasons discussed earlier. 



V. CONCLUSION 

We have performed molecular dynamics simulations for 
Xc layers on Ag(lOO) and Ag(lll) in order to determine 
the net sliding friction in the presence of a constant lat- 
eral external force. We find that fj depends sensitively 
on the corrugation amplitude uq of the Xe/Ag interac- 
tion potential and on the single-adatom microscopic fric- 
tion parameter ryy. For coverages approaching one mono- 
layer, the phonon- induced contribution r]ph{0) decreases 
rapidly, in agreement with the strong positive slope of 
the observed slip time r. A detailed study of riph{0) as a 
function of corrugation amplitude Uq suggests that, close 
to monolayer coverage, it is not possible to consistently 
describe both the magnitude and slope of the measured 
fj = 1/t without including a constant lateral electronic 
friction rjei ■ Our simulations indicate that this electronic 
contribution should be about 0.5 ns""'^ which is in qual- 
itative agreement with surface resistivity measurements 
and independent theoretical estimates. 

The overall picture that seems to be consistent with the 
available information is that the rapidly varying part of 
the observed sliding friction is due to excitation of atomic 
vibrations within the overlayer and that single-site elec- 
tronic friction processes contribute a roughly coverage- 
independent background. For a fully compressed Xe 
layer, phonon dissipation is weaker than the electronic 
friction. As the coverage is reduced, however, phonon 
processes rapidly become important and eventually, near 
the minimum of the observed slip time, dominate over 
the electronic dissipation channel. 

In summary, even though both experimental and theo- 
retical results on sliding friction of Xe on Ag have inher- 
ent uncertainties, our analysis indicates that both mech- 
anisms — phononic and electronic dissipation — do con- 
tribute. The relative weight of these channels depends 
strongly on coverage in the crucial region near a full 
monolayer. 
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